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Abstract 

We report a measurement of the ratio of branching fractions of the decays B ± — > J/ipiv^ and 
using the CDF II detector at the Fermilab Tevatron Collider. The signal from the 
Cabbibo-suppressed B^ 1 — ► J/ipir^ decay is separated from B ± — > J/ipK^ using the -► J/V^ 
invariant mass distribution and the kinematical differences of the hadron track in the two decay 
modes. From a sample of 220 pb~ x of pp collisions at y/s = 1.96 TeV, we observe 91±15 B^ — > 
J/i/m^ events together with 1883±34 B^ — > J/ipK^ 1 events. The ratio of branching fractions is 
found to be -► J/^jtt ± )/B(B ± -► J/ipK*) = (4.86 ± 0.82(stat.) ± 0.15(syst.))%. 

PACS numbers: 13.25.Hw 14.40.Nd 
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The B ± — > J I tpn^ decay is a Cabibbo-suppressed mode proceeding via a6-> ccd transi- 
tion. If the leading-order tree diagram is the dominant contribution, its branching fraction is 
expected to be ~ 5% of that of the Cabibbo-favored mode B ± — > J /ipK ± . Detailed predic- 
tions of the ratio are obtained using the hypothesis of factorization of the hadronic matrix 
elements , a theoretical approach widely used in the treatment of non-leptonic decays 
of B mesons. However, the absence of strong theoretical arguments supporting factorization 
and the use of phenomenological models, which are a source of theoretical uncertainties, 
weakens the reliability of those predictions, which need to be accurately tested on data. Un- 
til now, the measurements on the i? ± — > J j ipir^ decay were performed by many experiments. 
The BABAR collaboration reported B{B ± J/ifnr ± )/B{B ± -> J/ifrK*) = (5.37 ± 0.45)% 
with 244 ±20 B ± -> J/ijnr ± events [3|. The Belle collaboration reported B{B ± -> J/tjrn*) = 
(3.8 ±0.6) x 10~ 5 |4J. A previous study of the B ± — ► J/tpir^ 1 decay was also performed by the 
CLEO collaboration [5]. The result of this analysis supersedes the previous CDF result (|. 

This paper presents a measurement of the ratio of branching fractions B(B ± — > 
J/ip7T ± )/B(B ± — > J/ipK^. We use a sample of fully reconstructed B ± — > J/ipK ± de- 
cays, where J/ip — > corresponding to an integrated luminosity of 220 pb -1 of pp 
collisions at y/s =1.96 TeV collected by the CDF II detector at Fermilab between February 
2002 and August 2003. 

The CDF II detector is a multipurpose detector with a central geometry and has 
a tracking system surrounded by calorimeters and muon detectors. The components of 
the detector most relevant to this analysis are described briefly here. Charged particle 
trajectories are reconstructed in the pseudorapidity range \r]\ < 1.0, where r] = — ln(tan|) 
and 9 is the polar angle measured from the beam line 8]. Trajectories are reconstructed from 

n n 

hits in the silicon microstrip detector (SVX II) [9j and the central outer tracker (COT) [10J 
which are immersed in a 1.4 T solenoidal magnetic field. The SVX II consists of five 
concentric layers made of double-sided silicon detectors with radii between 2.5 and 10.6 cm, 
each providing a position measurement with 15 /im resolution in the r-0 plane. The COT 
is an open-cell drift chamber with 96 measurement layers, between 40 and 137 cm in radius, 
organized into eight alternating axial and ±2° stereo superlayers. The transverse momentum 
{Pt) resolution is u Pt /pt — 0.15%pr(GeV/c) -1 . Muon detectors consisting of multi-layer 
drift chambers are located radially around the outside of the calorimeter [ll|. The central 
muon detector (CMU) covers a range in pseudorapidity of \rj\ < 0.6. The central muon 



8 



extension (CMX) extends the pseudorapidity coverage to 0.6 < \r)\ < 1.0. 

The data sample used in this analysis required a dimuon trigger sensitive to J/tp — > [i + fi~. 
The CDF II detector employs a three-level trigger system to select events of interest effi- 
ciently At the first trigger level, muon candidates are identified by matching track segments 
in the CMU and CMX to coarsely reconstructed COT tracks obtained with the extremely 
fast tracker (XFT) [12]. Dimuon triggers use combinations of CMU-CMU and CMU-CMX 
muons with pt >1.5 (2.0) GeV/c for CMU (CMX) muons. For the data presented here, no 
additional requirements are made at the second level. At the third trigger level, a detailed 
reconstruction is performed and oppositely charged dimuon events with an invariant mass 
in the range of 2.7 — 4.0 GeV/c 2 are selected. 

In this analysis, we reconstruct — > J/ipK^ decays. B meson decay modes involving 
the well-known J/ip — > fi + n~ decay have been extensively used in other measurements at 
CDF, and their selection criteria are well established. We follow the selection requirements 
developed in the b hadron mass measurement 13J and apply them to the B decay mode 
of interest. To ensure the best momentum scale calibration, the data sample used for this 
analysis is also kept the same as that for the mass measurement. 

The B — > J / ipK^ 1 reconstruction begins by selecting J/ip — > candidates with pairs 
of oppositely charged tracks which satisfy the requirements of the dimuon triggers. J/if> 
candidates are further selected by requiring their invariant mass to be within 80 MeV/c 2 of 
the world average J/ ip mass [yj] . After a J/ip candidate is identified, any other charged track 
is assumed to be a kaon and is combined with the J/ip candidate to make a B^ candidate. 
The tracks of the kaon and two muons are then fitted to a common three dimensional vertex 
(3-D) while constraining the invariant mass of two muons to the world average J / ip mass [jjj . 
To ensure good vertex resolution, each track must have hits in at least three silicon vertex 
detector layers in the r — <fi plane and the probability resulting from the 3-D vertex fit is 
required to be greater than 1%. 

A number of further requirements are made to improve the signal-to-background sepa- 
ration. Prompt background, with tracks coming directly from the primary vertex, can be 
reduced by exploiting variables sensitive to the long lifetime of the B^ meson. To reduce 
prompt background, the transverse decay length (L xy ) of the B ± is required to exceed 200 
fim, where L xy is defined as the vector from the primary vertex to the B^ decay vertex 
projected onto the pr of the B ± candidate. To further reduce combinatorial background, 
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we require pt > 6.5 GeV/c for the candidate and Pt > 2.0 GeV/c for the hadron from 
the B^ decay. The values used in the above selection criteria are determined by an iterative 



optimization procedure in which the significance S/ V S + B is maximized. The quantity S 
represents the number of accepted signal events, in this case taken from a Monte Carlo sim- 
ulation sample, and B is the number of selected B^ candidates within the mass sidebands 
of the data. 

We measure the following ratio: 

B(B ± -» J/ip-K^) N J/i)1T ± e J/i)K ± 
B(B± -> J/$K±) N jm ± 

1 

= r obs x , (1) 

trel 

where r obs (= Nj/^±/Nj/^ k ±) is the ratio of the yields of each decay mode, and e re /(= 
e J/i>Tr±/ e J/i>K±) i s the relative reconstruction efficiency In this analysis, the quantity r obs is 
extracted from an unbinned maximum likelihood fit using the differences between the two 
decay modes in the mass distribution, and is corrected with e re i obtained from Monte Carlo 
simulation. 

To build the probability density function (PDF) used in the unbinned maximum likelihood 
fit, we choose the invariant mass of J/ip and a kaon (Mj/^k) as an observable. There are 
three components in the distribution of the Mj/^k variable: the B^ — ► J/ipK ± signal, 
the B^ — > J/ip-K signal and the combinatorial background. As demonstrated in the high 
statistics D and B mass reconstuctions with similar decay topology, the invariant mass 



1 described by a Gaussian function with 
[if] ]. Therefore, we model the 



3 Q W 



distribution of B^ — > J/ipK decay at CDF is we 
width determined by CDF's tracking resolution 
B^ — > J/i/jK^ 1 signal as a Gaussian centered at the mass of B^ 1 (Mb) with a width <jk- 
If the pion mass were assigned to the hadron track originating from the B ± — > J/ip^ 
decay, the resulting spectrum would be also a Gaussian centered at Mb- However, assigning 
the kaon mass to this track produces a spectrum partially overlapping the B^ — > J/ipK 
and shifted in the positive direction. The shifted invariant mass of B ± —>■ J/ib^ can be 

n 

calculated by an approximation, which has a good agreement with the exact value [17[ , 

M'l(a) ~ Ml + (1 + a){Ml - M%), (2) 

where Mx and M n are respectively the kaon and the pion masses. The purely kinematic 
variable a is defined as a = Ej/^ / P K , where Ej/$ is the J / ip energy and Pk is the magnitude 
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of the momentum of the hadron track. Using Eq. (j2j), the B — > J/tp7r signal is modeled 
as a Gaussian centered at Ai' B (a) with a width a^. We find ok and a n have almost the 
same value from the Monte Carlo simulation, so we constrain them to be the same value 
in the fit. We assume the background mass distribution is a first order polynomial. In the 
likelihood, we also include the PDF functions of a for B ± — > J/ipK and — > J/ipir^ 1 as 
the distributions for the two signals are found to be slightly different. We parameterize a 
PDF distributions from Monte Carlo simulation. We also parameterize the a distribution 
of the background, which is obtained from the mass sidebands of the data. These mass 
sidebands are chosen from 5.2 < Mj/^k < 5.24 and 5.4 < Mj/^k < 5.6 GeV/c 2 to avoid 
signal contaminations and other backgrounds from partially reconstructed B mesons that 
fall below 5.2 GeV/c 2 . The empirical functions used in the parameterizations are 

3 

hj/fx(at; fa, X i} a) = ^ fi( a ~ a)e" A ' Q (3) 
i=i 

3 

h bkg (a; fi, Aj, a) = f t (a - a)V Aia , (4) 

i=l 

where the symbol X denotes K or tc in Eq. ([3]) and f\, _/*2 and /3 are to be the fractional 
contributions of each type of function when the functions are properly normalized to 1. 
Because of the requirement on the px of the hadron track and also of the dimuon triggers, 
all a distributions show a cutoff around 0.5 in the a variable and these cutoff values are 
parameterized by a in Eq. ([3]) and (j4j). These parameters of the functions describing the a 
distributions are fixed in the fit. The a distributions of the signal and background, and the 
results of the parameters are shown in Fig. [TJ With models for each signal and background, 
and with the chosen observables, the PDF of the i th event is written as 



Pi =f t 



' G {M) m - M B , a) h jm {d) + ~^-G {M) m - M' B (a l ),a) hj/^a*) 



1 + r obs 
Hl-fs)B(M^ K )h bkg (a l ) 



(5) 



where f s is the fraction of signal events in the data sample and r a b s is the ratio between the 
yields of each signal. The functions, G (m i j/^ k — M B , crj and G (m)/^ k — •M.' B (a l ), aj , are 
Gaussians with a width a describing the mass distributions of B ± J/ipK^ 1 and B^ — > 
J/ifj7i ± , respectively, and B(Mj^ K ) is a first order polynomial function which describes the 
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FIG. 1: The a distributions of B^ — > J/ipK^, which are obtained with Monte Carlo simulation, 
and background obtained from the non-signal data sample. The solid curves are the corresponding 
parameterization functions from Eq. ([3]) and @. The a distributions of the two signals are very 
similar in shape due to the almost identical kinematics of the two decay modes. To avoid confusion 
from it, we plot the a distribution of — > J/i^K^ only. 

background mass distribution. The fitting range (5.2 < Mj/^ K < 5.6 GeV/c 2 ) is selected 
to avoid the backgrounds from partially reconstructed B mesons, but to include enough of 
the background region to determine accurately the background shape. C = YliLi Pi * s then 
maximized to obtain the best fit values for Mb, cr, f s and r & s . The fitter is extensively 
tested with Monte Carlo samples. 

The fit to 2683 candidates falling in the fitting range returns the signal fraction, f s = 
0.736 ± 0.012, and the ratio of the yields of each decay mode, r obs = (4.82 ± 0.81)%. These 
values give 1883±34 signal events in the B ± — > J/i/jK^ 1 decay mode and 91±15 events in 
the — > J/^ _± decay mode. The distributions in Mj/^k and a for the events in the data 
sample are shown in Fig. [2] and Fig. [3J along with the likelihood fit results. 

Possible biases in the fitting procedure are investigated by performing the fit on Monte 
Carlo samples generated by the PDF in Eq. (jSJ), with known composition and with the same 
size as the data sample. The difference of the ratio between the extracted and the input 
values is consistent with zero and the width of the pull distributions is one. 

In order to determine the ratio of branching fractions, the ratio of the yields of each de- 
cay mode must be corrected with the relative reconstruction efficiency. The relative recon- 
struction efficiency depends in turn on the different decay-in-flights and nuclear interaction 
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FIG. 2: The invariant mass distribution in the data sample (points) projected with the results of 
the likelihood fit; overall (solid line), B ± — > J/%l)K ± (dotted line), B^ — > J/^7r ± (dashed line) 
and background (dashed-dotted line). The inset shows the magnified region of the B ± -► J/V'vr ± 
signal. 




FIG. 3: The a distribution in the data sample (points) compared with the results of the likelihood 
fit; overall (solid line), B ± — > J '/ipK ± (dotted line), B^ — > J/ipir^ 1 (dashed line) and background 
(dashed-dotted line). 



probabilities of the kaon and pion from the two decay modes and on the slightly different 
track momentum spectra. The relative reconstruction efficiency for the two decay modes is 
e re i = 0.991 ± 0.005 which is derived from the Monte Carlo simulation. 

In this analysis, we use a Monte Carlo simulation to parameterize the a distributions of 
each signal and to determine the relative reconstruction efficiency for the two decay modes. 
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The Monte Carlo generation proceeds as follows. Transverse momentum and rapidity distri- 
butions of single b quarks are generated based on next-to-leading order (NLO) pertubative 



QCD calculation 
fragmentation 



181 ] . B meson kinematic distributions are obtained by simulating Peterson 



191 ] on quark-level distributions. Additional fragmentation particles, corre- 
lated bb production and the underlying event structure are not generated. The B meson 
spectrum used in the Monte Carlo simulation is from the inclusive B — > J/ibX measure- 

n n 

ment [Tj. The cleomc program [20] is used to decay B ± mesons into the final states of 
interest. The simulation of the CDF II detector and trigger is based on a geant 2l| 
description. 

Since both decay modes of interest have almost identical decay topology and kinematics, 
most systematic uncertainties cancel in this ratio measurement, including uncertainties in 
total integrated luminosity and trigger and reconstruction efficiencies. Remaining systematic 
uncertainties come from the uncertainties in the shapes of the mass distribution, the param- 
eterized PDFs in the a variable, and from the determination of the relative reconstruction 
efficiency. The largest systematic uncertainty originates from the unknown shape of the 
combinatorial background in the mass distribution. To estimate this effect, a second order 
polynomial function is considered as an alternative model for the shape of the background 
mass distribution. The modelling of the width of the invariant mass distribution is deter- 
nnned fr0 ,„ —ntn scale reso.ntion stndies Q An alternative ,„ode, tan a staple 
Gaussian is to include an additional Gaussian for potential different momentum resolutions 
of tracks reconstructed in different detector geometry coverage. We replace a Gaussian with 
a double Gaussian for modeling each signal mass distribution and fit again to evaluate the 
uncertainty coming from the non-Gaussian tails in the B^ — > J/ipK ± mass distribution. 
The uncertainties in the function parameters describing the a PDFs , in Eq. ([3]) and (jlj), 
generate an uncertainty for the ratio measurement. The contribution of this uncertainty is 
estimated by performing the fit by varying the parameterzation variables of the PDFs by 
the ±1(7, obtained from Monte Carlo simulation. The uncertainty in e re / originates from the 
uncertainties of the nuclear interaction and the material description in the detector simu- 
lation. The geant simulation calculates nuclear interaction probabilities of ~ 4% for 7r + , 
7r~ and K~ , and ~ 3% for K + . We then assign a 25% uncertainty to the calculated nuclear 
interaction probabilities as the uncertainty of the detector material description in the de- 
tector simulation, and take the resulting uncertainty in e re i as a systematic uncertainty. We 
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determine the total systematic uncertainty of 3.0% on the measurement by adding the indi- 
vidual uncertainties in quadrature, and the contributions from each source are summarized 
in Table [H 



TABLE I: Summary of systematic uncertainties for the ratio of branching fractions, B(B 
J/tp7r ± )/B(B ± -» J/ipK^). 



Source 


Uncertainty of 




the ratio (%) 


Background Shape 


2.5 


Non-Gaussian tail of — ► J/ipK 


1.2 


a PDFs Parametrization 


1.0 


Relative Reconstruction Efficiency 


0.5 


Total Uncertainty 


3.0 



From Eq. flTJ), we derived the ratio of branching fractions, 

B(B±Zj/Tk±) = (4 ' 86 ± °- 82 ( stat ') ± 0-15(syst.))%, 

where the first error is statistical and the second is systematic. 

In conclusion, we present the measurement of the ratio of branching fractions between 
5 ± — > J/tpTi^ and —>■ J/tpK ± . This result is consistent with theoretical expectations 
and the previous measurements, and will improve the present world average (5.3±0.4)% jl^j ]. 
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